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Abstract: 

We employ nonperturbative flow equations for an investigation of the effective action 
in Yang-Mills theories. We compute the effective action T[B] for constant color mag- 
netic fields B and examine Savvidy's conjecture of an unstable perturbative vacuum. 
Our results indicate that the absolute minimum of T[B] occurs for B = 0. Gluon 
condensation is described by a nonvanishing expectation value of the regularized 
composite operator F^^^F'^'^ which agrees with phenomenological estimates. 



1 Introduction 



A perturbative calculation of the effective action for constant color-magnetic fields 
B in Yang-Mills theories indicates that the configuration of lowest Euclidean action 
does not correspond to vanishing fields This has led to many interesting 

speculations about the nature of the QCD vacuum. Unfortunately, perturbation 
theory is clearly invalid in the interesting region in field space. It breaks down both 
for vanishing magnetic fields and for fields B corresponding to the minimum of the 
perturbatively calculated effective action T[B]. A nonperturbative computation of 
T[B] was lacking so far. In this paper we make an attempt to evaluate T[B] based on 
the non-perturbative method of the average action 0. Within our approximations 
we find that the minimum of T[B] occurs for a vanishing magnetic field i? = 0, in 
contrast to the perturbative result , • This strongly suggests that the minimum 
of the effective action occurs for vanishing gauge fields = 0. This implies a 
vanishing expectation value < >= and Lorentz invariance of the QCD vacuum 
is trivially realized^. 

The vanishing expectation value < A^ > does by no means indicate that the 
QCD vacuum is simple. Gluon condensation phenomena may be described by non- 
vanishing expectation values of composite operators such as < F^,,F^^ >. We also 
discuss this issue in the context of the average action and compute the expectation 
value of a suitably regularized operator. Our result for the gluon condensate agrees 
with phenomenological estimates. (The anomalous trace of the energy momentum 
tensor induced by the condensate < F^yF^" > is found to be —0.015 GeV*). This 
may be considered as a first analytical computation of this condensate from "first 
principles". However, in view of the large possible truncation errors and the un- 
"'^See ref. for a discussion of possible realizations of Lorentz symmetry for < >^ 0. 
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certainty in the precise correspondence between the computed operator and more 
phenomenological concepts of a gluon condensate it seems premature to draw any 
quantitative conclusions. We rather consider this work as an encouraging first step 
towards a quantitatively more reliable computation of the gluon condensate. 

The average action is the effective Euclidean action for averages of fields which 
obtains by integrating out all quantum fluctuations with (generalized) momenta 

> A;2. It can be viewed as the standard effective action T computed with an 
additional infrared cutoff ~ k for all fluctuations. In the limit — > the average 
action equals the usual effective action, Fq = F. For A; > no infrared divergences 
should appear in the computations and a lowering of k allows to explore the long- 
distance physics step by step. The fc-dependence of the average action is described 
by an exact nonperturbative evolution equation The structure of this equation 
is close to a perturbative one-loop equation but it involves the full propagator and 
vertices instead of the classical ones. For gauge theories it can be formulated in a 
way such that Ffc[A] is a gauge-invariant functional of the gauge field A 0, 0, 

The flow equation for the pure non-abelian Yang-Mills theory has already been 
solved with a very simple approximation - the average action F^ has been trun- 
cated with a minimal kinetic term ~ Z^F^j^F^" . From the /c-dependence of 
the running of the renormalized gauge coupling g{k) has been derived for arbitrary 
dimension d. Most strikingly, this lowest-order estimate suggests that the wave 
function renormalization Z^ reaches zero for k = k^ and turns negative for k < k^. 
Here the scale k^o can be identified with the confinement scale, i.e., the scale where 
the renormalized gauge coupling g'^{k) ~ Z^^ diverges. This was interpreted as an 
indication for the instability of the perturbative vacuum with = 0, similar to 

^For the relation to earlier versions of exact renormalization group equations |^ see ref. . 
•^See ref. jl^, for alternative formulations. 
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the perturbative result IQ. Indeed, negative would imply that the minimum 
of Ffc must occur for nonvanishing gauge field for all k < [|12]. It is obvious, 
however, that a truncated effective action ~ Z^B^ is insufficient to describe such 
phenomena. One expects positive Tk for large B and this cannot be accommodated 
with a negative in the "lowest order truncation" . 

In this paper we enlarge the "space of actions" by considering for an arbi- 
trary function of constant magnetic fields, i.e., Tk ~ Wk{\B'^). The flow equation 
describes how the function Wk changes its shape, starting from a linear dependence 
Wk = \ZkB'^ for large k. The approximations (truncations) employed for the com- 
putation and solution of the flow equation are summarized at the end of sect. 2. 
In particular, we solve numerically the flow equations in a truncation where Wk is 
approximated by a polynomial of order B^. We find that in contrast to the pertur- 
bative computation of ref. the minimum of Wk remains at i? = for all k. The 
vanishing of Zk for k^o > turns out to be an artifact of the lowest order truncation. 
In the truncation of order B^ the gauge coupling grows large but remains finite. In 
this approximation the ground state value of the gauge field vanishes, in contrast 



to speculations in ref. [0. 

For the description of gluon condensation in terms of the expectation value of a 
suitably regularized gauge-invariant composite operator F^^F^" we use exact flow 
equations including composite operators [|l3l. Following this approach we enlarge 
in sect. 4 our truncation by studying a general function Vk{\B'^,x), where x is a 
composite field related to the operator F^yF^^^ . We find that x develops indeed an 
expectation value which is stable for 0. This implies a nonvanishing gluon 

condensate. 

In the following we shall consider the flow equation for the pure SU{N) Yang- 
Mills theory. The general formalism of our formulation with ghost fields and back- 
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ground fields A is briefly reviewed in appendix A, to which the reader could turn at 
this point. There we write down the corresponding modified Slavnov- Taylor identi- 
ties fill] and a further identity which describes the dependence of the effective action 
on the background field [jl4|, [|15]. We also give in this appendix the full exact 



flow equation including the ghost sector. Since we are interested in the ground state 
and there exist always classical solutions with vanishing ghost fields, we concentrate 
here on the zero ghost sector. Then the effective action reduces to a functional 
rfc[y4, y4] depending only on the gauge field A and the background field A. It is 
invariant under simultaneous gauge transformations of both fields. With a suitable 
truncation for the ghost dependence of the full effective action (cf. appendix A) 
one obtains the following renormalization group equation for the scale dependence 
of rfe[A,A] {t = lnk): 



-Tr„ 



Tf[A, A] + Ru{Tf[A, A])y' ^R,{Tt\A, A]) 
-D^[A]D,[A] + R,{-D'[A]))-^^R,{-D'[A]) 



Here the first trace on the r.h.s. arises from the fluctuations of the gauge field. It 
involves on integration over space-time ( "x" ) as well as a summation over color ( "c" ) 
and Lorentz ("i^") indices. The "color" trace is in the adjoint representation. The 
second trace is due to the Faddeev-Popov ghosts and is over space-time and color 
indices only. The precise form of the infrared cutoff is described by the function 
Rk- It is convenient to choose a smooth cutoff which vanishes for large covariant 
momenta 

Rk{u) = u[exp{Z^\/k'^) - 1]^^ (1.2) 

For the purpose of comparison we will also consider in this paper a second choice 
for the cutoff, namely simply a constant 

Rk = Zkk'' (1.3) 



In eqs. ( |1.3| ) the wave function renormalization Z]^ could be a matrix in the 

space of fields which may even depend on A. (In ref. we used 2^ = 1 for the 
ghosts and a fc-dependent constant Z]. = Zp^k for all modes of the gauge field.) 
Eq. ( p. ■ 1|) is a functional differential equation where r^^^[A,A] denotes the matrix 
of second functional derivatives of with respect to A, with the background field 
A kept fixed. The modes of the gauge field are declared "high-frequency modes" 
or "low-frequency modes" depending on whether their eigenvalues with respect to 
the operator [A, A] = vf [A, A] J are larger or smaller than fc^, respectively. 
Therefore it is this operator which appears in the argument of and in this sense 
Rk acts as an effective infrared cutoff by suppressing the low frequency modes. 

In this paper we aim for an approximate solution of the flow equation by a 
truncation of the most general form of Tk[A, A\. We derive in sect. 2 the evolution 
equation for the effective action for constant color magnetic fields B. This equation 
is further approximated in sect. 3 by using a polynomial ansatz for Tk{B) of order 

and solving it numerically. In sect. 4 we enlarge the truncation by introducing 
a field for the composite operator F^^F^^ . This will permit the investigation of the 
flow of the expectation value of F^j^F^^ and an extraction of the gluon condensate 
for A; ^ 0. Sect. 5 finally contains our conclusions and a brief discussion of further 

possible developments. Since some of the necessary computations are technically 

^In we used the classical 5'^^ [A, A\ rather than F^^"* [A, A\ for this purpose. While preserving 
all the general properties of Tk, the new flow equation is much easier to handle from a technical 



point of view. This change also entails the different positioning of Zk in (1.2) relative to the one in 



for the simple truncation used there one has r^^^ = Zp^kS'^^^ for part of the modes. We only 



will use the definition (1.2) if for high momenta cp' ^ oo one has u — > oo. We should mention, 
however, that for negative eigenvalues of ' [A, A] the vanishing of Rk for A; is guaranteed only 



if the ratio u/k^ in (1.2) remains finite. This problem concerns mainly the approach to convexity 
of the effective action for /c — > and is of no relevance for the present work. 
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involved, we display a lot of this material in various appendices. 



2 Evolution equation of the effective action for 
constant color- magnetic fields 

In order to find nonperturbative approximative solutions of ( |1.1| ) we employ the 
following ansatz for F^: 

T,[A,A] = J d'x W,i\F^^ix)Frix)) + d'xj:iD,[A]iA^ - An)l (2.1) 

Here Wk is an arbitrary function of the invariant jF^ with F^^^ the field strength of 
the gauge field A. The fc-dependence of Wk will be determined by inserting (|2.1|) 
into the evolution equation ( |1 . 1| ) . If we think of Wk{9), 

e ^ \F^^Fr, (2.2) 

as a power series in 6 our ansatz contains at this point still invariants of arbitrarily 
high canonical dimension. Since all invariants which occur are of the form 6\ only 
the dimensions = 1,2,..., actually occur in this truncation. Also, for a fixed 
dimension 4/ the truncation (p.l|) does not contain a complete basis of operators. 



Nevertheless one may hope that (|2.1|) gives a qualitatively correct picture of the 
effective action for constant color magnetic fields in the regime where the renor- 
malization group evolution has already drastically modified the classical Lagrangian 



1 f2 



F . We remark that effective actions which depend on 9 only play also a central 



role in the leading-log models [|T7| of QCD. 

The second term on the r.h.s. of ( |2.1| ) is a standard background gauge- fixing 
term ^ with a fc-dependent gauge- fixing parameter a^. In addition to the 
/c-dependence of the function Wk{0) we should, in principle, also compute the k- 
dependence of a^- We will omit this here since the general identities of appendix A 

6 



imply that a is independent of in a first approximation. Furthermore, we will see 
that within the truncation and for a suitable choice of the infrared cutoff 
the evolution equation for Wk becomes independent of a^. 

We should mention at this place that a truncation is actually not completely 
defined by the terms retained but rather by specifying which invariants in the most 
general form of are omitted. One may parametrize a general F^ by infinitely 
many couplings multiplying the infinitely many possible invariants which can be 
formed from the gauge fields consistent with the symmetries. In the corresponding 
infinite dimensional space a truncation is a projection on a subspace which is defined 
by setting all but the specified couplings to zero. (In our case the subspace remains 
infinite dimensional.) In practice, we will choose a particular test configuration 
corresponding to a constant magnetic field. The truncation should then be under- 
stood in the sense that we use a basis for the invariants where all invariants except 
those used in ( |2.1|) vanish for the test configuration. By putting the coefficients 
of all invariants which vanish for the test configuration to zero the truncation is 
uniquely defined. A computation of therefore amounts to a computation of the 
fc-dependent effective action for a (particular) constant magnetic field. 

As a particularly convenient test field we choose a covariantly constant color- 
magnetic field 1^ with a vector potential of the form 

A^(x) = n^A,(x) (2.3) 

Here ri^ is a constant unit vector in color space (n^n^ = 1), and A^(x) is any 
"abehan" gauge field whose field strength 

F^v = d^,K - d^k^ = Be^^ = const (2.4) 

corresponds to a constant magnetic field B along the 3-direction, say. (We de- 
fine = —^21 — with all other components vanishing.) Hence we have 6 = 
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In summary, the solution of the evolution equation will produce the effective 



action for constant color magnetic fields of the type ( |2.3| , |2.4|) , where an additional 
infrared cutoff k is present. One is finally interested in the limit A; ^ 0. For the 
computation of this effective action we make approximations which amount to the 
following truncations of the effective action: 

1. The ghost sector is approximated by its classical form, as discussed in appendix 
A. 

2. The remaining gauge field dependence of is approximated by 

3. The resulting flow equation which will be derived in this section is a nonlinear 
partial differential equation for a function of two variables, Wk{9) = W{k,9). 
We solve numerically in sect. 3 only an approximate version of this differential 
equation where Wk{0) is truncated to a polynomial in 6 of degree three. 

We do not expect these approximations to yield a quantitatively precise result 
in the range of k where the gauge coupling is large. Nevertheless, it remains an 
interesting question if the speculated effect of an unstable perturbative vacuum 
persists in this picture. Our findings indicate, in contrast to perturbation theory 
[§] , that the minimum of the effective action occurs for = 0. We emphasize that 
despite the approximations made our approach goes far beyond the perturbative 
calculation of refs. We conclude that there is no reason to believe that the 

configuration B = {or = 0) is unstable. 

The evolution equation for Wk{6) is computed for arbitrary dimension d in ap- 
pendix B. One finds 



+Vd 



\ --1 



Wl + B^Wll Wl) \W'J 
-V, rdxxi-\H (^)-h (-\ 

Jo \CikJ \CikJ 

with = {dWk/de){lB^), etc., = J d'^x, and 



/ dxx^~^H{x) 
Jo 



-1 ^ 2'^+V2r 



(2.5) 



(2.6) 



Here we have introduced the convenient abbreviation (for = 



H{u) = {u + Rk{u))-'—Rk{u) 



(2 + |ln^.) 



exp 



{2 + fjnZk)Zke[u + Zkk 



21-1 



for (1.2) 



for (1.3) 



(2.7) 



In the second, ghost-type, trace the function Hg{u) is defined similarly, but with a 
different factor 2^ = 1. The last trace accounts for a possible difference between Z^ 
and Zk in ( |B.20| ). Here H is obtained from H by replacing Zk by Z^- 
The eigenvalues of the operator 

Vt = -D^ + 2igF 

are known explicitly 0. They are parametrized by a (c? — 2)-dimensional momentum 
qf^ which "lives" in the space orthogonal to the 1-2 plane, and a discrete quantum 
number n = 0,1,2,... which labels the Landau levels. The spectral sum for the 
function H{x) = H{Wl,x) reads^ 



271 ^„J (27r)* 



^The momentum integration is absent for d = 2. 



{d-2)H{q^ + {2n + l)g\ui\B) 
+H{q^ + i2n + 3)g\ui\B) 

+Hiq^ + {2n-l)g\ui\B)} (2.8) 

Here z/;,/ = 1,...,A^^ — 1 are the eigenvalues of the matrix n^Tz in the adjoint 
representation. We note that for n = and sufficiently small the eigenvalue 
— gliyilB in the third term on the r.h.s. of (|2.8|) can become negative. This 



instability [y, |I8[ causes severe problems if one tries to compute the standard one- 
loop effective action in the background of a covariantly constant magnetic field [l^ . 
In our approach this problem is cured by the presence of an IR regulator. Eq. (|2.8|) 
can be rewritten {d > 2) as 



Q''Tt,,l[H{Vt)] = ^ E 9Wi\B r dxx 

TT Jo 



dJ2H{x + {2n + l)g\vi\B) + H{x - -g\vi\B) - H{x + g\iyi\B) (2.9) 



n=0 



but it cannot be simplified any further in closed form. The other traces in ( |2.5D are 
given by 

„. N^-l oo 

Q-'Tt,,[H{-D^)] = ^ E 9Wi\B / dxx-^~^ E H{x + {2n + l)g\i^i\B) (2.10) 

^ 1=1 -^0 n=0 

and similarly for H and Hq- 

In this paper we use two different methods in order to (approximately) com- 
pute the spectral sums (p.9|) and (|2.1CI|) . In appendix F we shall represent them as 
Schwinger proper-time integrals [0. This method leads to compact integral rep- 
resentations which are valid for all values of B, but it has the disadvantage that it 
works only for the cutoff function Rk{x) defined by ( |A.31| ) which leads, as we shall 
see, to ultraviolet problems. 

The second method consists of expanding the r.h.s. of (|2.9| ) in powers of B. 
It is applicable if gB <^ . It works for any function Rk{x) such that ultraviolet 
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problems can easily be avoided. The condition gB <^ k"^ guarantees that we may 
express the sum over n by an Euler-McLaurin series, and that the terms become 
small rapidly. We will concentrate mainly on the second approach. In this manner 
rof) turns into 

= {N' - l)^[WUlB')]-i f dx fdyxi-'H{x + y) (2.11) 



^ m=l \ 1=1 



with 



Here B2m are the Bernoulli numbers. In a second step one has to expand the 
-B^-dependence of W^,. Note that only even powers of B occur in this expansion. 
The group-theoretical factors Y^f^i"^ are discussed in appendices C and D. In 
particular for SU (2) these factors equal 2 for all values of m. 

If we use the Euler-McLaurin series ( p.ll|) and a similar expansion for the trace 
( CT) in eq. (gj), we find for SU{N): 



2n 



-2m 



TT 



m=l 

d 



with Tm defined in appendix D and the constants Em given by 

i^m = ^(2^"^-^ - l)52m. (2.14) 
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The dimensionless integrals 



d,m 



{Zkk 



X) 







/ , \ 2 m— 1 
CtX X2 — 



dx 



d 

ax 

— 1 



(^fcA;2)-2 / dxx^-^Hix) 

2 — 7] Jo 

2 ^ /"OO /-oo ^ 

{Zkk^)-2 dx dyx2-^H{x + y) 

2 — 1] Jo Jo 

r dx r dy ^^"'(^ + ^^ 
JO Jo exp(x + y) — 1 



(2.15) 



occur as a consequence of ( |2.7|) . The second equahty in eqs. (|2.15|) uses the expo- 
nential cutoff (|1.2| ). For this choice^ we note for later use that in 4 dimensions 



4,m, 



B2m-2, r\ = 4 = 2C(3) 



(2.16) 



where C, denotes the Riemann zeta function. The evolution equation ( p.l3|) is the 
central result of this section. It constitutes a partial differential equation for a 
function of two variables, W{9, k). 

In eq. ( p.l3|) we have introduced the anomalous dimensions 



In Zfc 

dt 



(2.17) 



A convenient choice for the wave function renormalization constants used in the 



^In contrast, the quantities rf and r| are not well defined for the choice (A. 31). The ultraviolet 
divergence indicates an incomplete "thinning out" of the high momentum degrees of freedom for 



a simple mass like infrared cutoff. Even though eq. ( ^.13 ) was derived by choosing a specific 
background field, this evolution equation does not depend on the background we used for the 
calculation. By employing derivative expansion techniques pl|, p3] it should also be possible to 



derive (2.13) without ever specifying a background. In the case at hand the method presented here 
is by far simpler, however. 



12 



infrared cutoff Rh is 



Wl{Q) for k>k 



np 



WL (0) for k<k 



np 



Zk = — (2.18) 

where knp is a typical momentum scale which characterizes the transition from the 
perturbative to the nonperturbative regime. We account for the possibility that 
Wl.{0) may turn negative for k smaller than a "confinement scale" k^o, whereas 
must always be strictly positive. More precisely, if \dlnWl.{0)/dt\ becomes of order 
unity for small scales k, we choose knp to be the scale where |?7(A;„p)| = 1.5. For 
k < knp all couplings run fast anyhow, and an improvement of the scaling properties 
of Rk by the introduction of a fc-dependent wave function renormalization seems 
not necessary. The choice Zk = guarantees that the infrared cutoff acts on the 
longitudinal modes in the same way as on the transversal modes. It implies that the 
flow equation for -^WkiO) becomes independent^ of ak and therefore independent 
of the "gauge fixing" in our truncation! In this paper we can therefore neglect the 
running^ of a^, and one has fj = 0. 

It is convenient to express the flow equation in terms of renormalized dimension- 
less quantities 

= k'-'Z^Y 
^ = g^k-'^ZkO 



0- independent constant in Wk is irrelevant. 
^This can be inferred from a first-order approximation to the solution of the general identities 

which govern the dependence of rfc[^,yl] on the background field A Q, jf^. 
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Wk{^) = g'k-'^Wkie) (2.19) 

Switching to a notation where dots denote derivatives with respect to instead of 
9 and with d/dt now taken at fixed -d one obtains 

d 

— = -(4 - ri)wk{'d) + Mwki^) 

+(2 - v)v,g'{M^))-'^ | (ci-iKd-2) ^^, _ 



+2(d - 2)^;rf52 ^ 7-^E^ro''"X2i?)'" + const (2.20) 

m=l 



This nonhnear partial differential equation for the function w{'&,t) does not show 
an explicit t- (or A;)-dependence of the right-hand side any more. 

One needs in addition the running of the renormalized gauge coupling g and the 
anomalous dimension rj, which are related by 

P,^^^^{d-A + v)g' (2.21) 

For k > knp we have by definition Wfe(O) — 1 and r) can be determined by 

^«;fc(0) = = r/ - 2(d - 2)v,4''ng'' ((2 - r])Ct - (4 - r])E^) 

-(2 - r,)v,g' (24 + ^(^^llM^^N' - 1)4^ *,(0) (2.22) 

With n ^ N,C^ ^ ^-2,Ei^ this yields 

- (^jVd{d-2){26-d)4''g^-2hdg^W2^ 

(1 - jva{d - 2) (25 - d)4'g' + ^^2) (2.23) 



where 



(2.24) 
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W2=Wk{^) (2.25) 

For general d the constants Tq'^, and depend on the precise choice of the infrared 
cutoff except for d = A where Tq'^ = 1 is cutoff independent. The running of the 
renormahzed gauge couphng g is now fully determined by eq. (|2.21| ) . It depends on 



the additional couphng W2 ( p.25| ) which will be discussed in more detail in the next 
section]^. 

Specifying the initial value g'^{A.) and the function w\{i)) at some high momentum 
scale A the form of Wk{'&) and g'^{k) are completely determined by the flow equation 
( |2.20| ). Solving for k the function Wo('(9) specifies the effective action in our 



truncation. If necessary, one has to replace for k < k^p eq. ( p.23| ) by ?7 = 0,g^{k < 

knp) 9 {knp)' 

Before closing this section, we briefly comment on the range of convergence of 
the Euler-McLaurin series in our case. For large m one has 



lim Ci=d lim E^a ~ vr"^™ (2.26) 

m— >oo m— >oo 



For iV = 2 and = 4 we find (cf. ( glBD ) that the coefficients of (2^9)*" in eq. (12^201) 

diverge ~ vr~^™i?2m-2 ~ 7r~^''"2~^™(2m — 4)!. For small nonvanishing d the first 

terms of the series have a very rapid apparent convergence, but the series finally 

diverges due to the factorial growth ~ (2m — 4)!. We can therefore safely use this 

series only for the derivatives w'^^\'d = 0) with finite n where convergence problems 

are absent since only a finite number of terms in the sum contributes. The situation 

is probably similar for N > 2 and /or d 7^ 4 as well as for many other choices of the 

infrared cutoff. In contrast, the original sums over n in eqs. (|2.9| ), ( p.lO|) always 

converge since for B > the contributions from sufficiently high values of n are 
^For W2 = we recover the result of rcf. |^] except for the factor (25 — d) in the denominator 



in ( 2.23 ) which was (24 — d) previously. This difference is due to a slightly different choice of the 



Z-factors in the infrared cutoff. 
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exponentially suppressed. An explicit evaluation of these sums is possible for the 
simplified masslike IR-cutoff (|A.31| ). This is described in appendix F where we will 
also see the reason for the ultraviolet divergence of rf and for this particular 
cutoff. 



3 Polynomial truncations 

In this section we concentrate on c? = 4 with the exponential infrared cutoff (1.2). 
One could solve the flow equation ( p.20| ) numerically. Instead, we further simplify 



here the truncation in order to get a first idea of the physical contents of ( ^.1^ ). We 
include in Wk{'&) only terms which are at most quadratic in -d: 

2 o 

Wj{k) = (^-^^ Wk\^=o (3.1) 

Thus the truncation for Tk is parametrized (up to an irrelevant constant) by three 
couplings: 

4 g\k) ' 32 g\k)k^ 

384 c/2(A;)A;S 

The short distance or "classical" theory is specified for A; = A by (Za = 1) 

g\t.) = -g\ w7i(A) = l, «;2(A) = «;3(A) = (3.3) 
For k > knp it follows from the definitions (|2.18 ) and (|2.19|) that Wi{k) = 1, and 



T,\A A] = fd'xl-^^^^F^ ]_ g'w2{k) . 

+ i^te^r)1 (3.2) 



g (k) is determined by eqs. ( |2.21| ), (|2.23| ). For k < knp we use g (k) = g {knp) 
instead and keep 

wiik)=WkiO) (3.4) 
16 



as the independent running coupling constant The flow equations for the partial 
derivatives Wj{k) = wl^\i!} = 0) follow by differentiating eq. (|2.2CI|) with respect to 
{}. For example, one has 

d 

-lb = rjw + Miv — (2 — vi)vdg^w~'2 ■ 



{d-2)J2 ^m{Ci - E^) 4''^{2dw^)'^-\2mw'' + (4m - d)dww) 



dt 

oo 

'''d rp \ „d,,m /r,„Q„-2\m-l / r,^„-,2 
k m=l 

^^^y w + 2§w {w + 2§wy w{w + 2^w) ^ 

Ad{d-l){d-2){N'-l)Ti- 

4 w 



+A{d - 2) W E mT^E^rt''^{2^r-^ (3.5) 

m=l 



We observe that Wq does not appear on the r.h.s. of the evolution equations for 
Wjij > 1, and we omit in the following this irrelevant constant. 

For k > knp the evolution equation for the running gauge coupling reads {d = 4) 



{UN - 3H^W2) 



dt 247r2 

whereas for k < knp one uses the flow equation for wi 



(3.6) 



d g\knp) (UN 
a7^i = Q 2 ~^ (3.7) 



where = IGn'^hi = rj + 3{N^ - 1)4 = 2{3N^ - 2)C(3) for the choice It 
is obvious that for W2 < the gauge coupling g^{k) always increases ( |3.6| ) until at 
k = knp the anomalous dimension \ri\ reaches 1.5. If W2 remains negative for k < knp, 
the coupling wi decreases until it reaches zero at the confinement scale fcoo > 0. (If 
we define g'^{k) = g'^{knp) /wi, this coupling diverges at the confinement scale.) The 
issue is different for W2 > 0: The coupling Wi does not reach zero since for small 
enough wi the second term in eq. (|3.7|) would cancel the first term. One therefore 
needs an estimate of W2{k). 
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Let us first consider the regime k > knp where the evolution equation for W2 
reads (with E2 = -j^, C*! = 

^ (A^ ^ 9' f 4,2/127 29 \ 

+ (2 - v) (Srt + liN'- 1)4^ wl -{2-v) (r^ + ^(iV^ - l)^^') ^3 

We observe the appearance of the couphng 1^3 whose evolution is given by (for details 
see appendix F) 

9 /ON 1 /442 137 \ 

-^3 = (8 + V)ws + - ^V) rs (3.9) 

+^(2 - r])T2W2 - 6C(3)(2 - vW2N^ + 23)w'l - {9N^ + 8)^2^3 + N^w^^ 

We neglect the term ~ W4 which is consistent with our approximation. In the 
perturbative region, where Ng'^/lGn'^ is small, it is easy to infer from ( p.8| ) that W2 
is of the order g"^. In fact, we may neglect rj and the terms ~ and ~ W3 in the 
curly bracket in (|3.8| ). In lowest order one finds for the ratio ^ an infrared stable 
fixed point: 

d ( W2\ W2 , 127 4,2 in 



which is approached very rapidly. (The last equality in ( p.ll|) holds for = 2 and 
uses Tg'^ = I for the exponential cutoff.) Similarly, one obtains the lowest order 
fixed point 

22lT,g\k) 

= 37800T6^ ^'-^'^ 
Actually, as we show in Appendix F, all ratios Wnj g^ reach perturbative fixed points 
for > 3. This justifies the approximation (|3.1|) at least for small enough g^{k\ 
(We observe Wn ~ ^2(n-i)^-4(n-i)_) 
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It is interesting to insert the value ( p. Ill ) into the /9-function for eq. ( ^3.61 ). 

Expanding in powers of g'^ one has 

dg^ 22iV / 77iV2^6 

-m = -^16^ - ^ (16^ + - 2)C(3)— 
22 g'N /rr 127 / 2 \\ 

We note that without the term ~ W2 the coefficient ~ g^ exceeds the perturbative 
two-loop coefficient oiily by a little more than 10%. It is recomforting 

to find the contribution from W2 in the same order of magnitude as this difference. 
We emphasize that a full computation of (3g2 in order g^ should take additional 
invariants into account, as for example [FF^ or {D^F^'^Y. We also observe that 
for k = knp, i.e., for r/ = — | one has approximately Ng"^ /IQv:'^ = | so that the 
validity of perturbation theory extends roughly to all k > knp. 

Let us finally consider the regime k < k^p- Our truncation permits, in principle, 
a first-order type transition where the absolute minimum of w jumps from = to 
a nonzero value > 0. For positive the polynomial ( p.l| ) is bounded from below 
for ?9 > 0. For wi > there is always a local minimum at ^9 = 0. Two additional 
extrema are present if 

wl > 2wiWs (3.14) 

There is a minimum for positive -^q 



-W2 + ^jwl-2wiw^ 

Vo = (3.15) 

and, for Wi > 0, a. maximum at 

-W2- ^Jwl-2wiW3 

Vmax = (3.16) 

For Wi < the origin = turns to a local maximum. The critical set of couplings 
where the two minima are of equal height {wk{'&o) = 0) corresponds to 

wi = -— 3.17 

8^3 
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For W3 smaller than the critical value, the absolute minimum occurs at •&(). This 
would correspond to the picture where the perturbative vacuum is unstable, i.e. 
< >^ 0. 

For k > knp where Wi = 1 an inspection of the flow equations (cf. ( p.l2|) ) 



shows that 2w3 remains larger than There is therefore only one minimum at 
■(9 = 0. If W2 remains negative for k < knp, the coupling wi would decrease towards 
zero. If furthermore W3 stays positive, the condition ( p.l7|) would then necessarily 



be met for some scale k > k^. On the other hand, if W2 turns positive, Wi will also 
remain strictly positive, and the only minimum occurs at i} = for all values of k. 
We have solved numerically the system of flow equations for k < knp 



d , , g^Kp) f /29 4 7 ,\ 



2 

+4C(3)[(9iV2 ^ _ (3^2 _ 

wf Wi ' 

d . , 9' (knp) f T-3 /137 6 31 2\ 87 

7/!o = S7/!n -\ ^ — ^ 7/ir -\ 1117 -\ 1 



— = + \-^[ tttf^i + ttt^i ) + y^'^2w?w2 

together with eq. ( p.7|) for Wi. (According to our truncation, we put W4 = 0.) 
For = 3 we find that remains positive for all fc, whereas W2 changes sign for 
k ^ 0.5A;„p ~ Aqcd- In consequence, all three couplings run to fixed values as /c — > 0, 
given by 

= 0.13, = 4■10-^ w;3* = 4-10-^ (3.19) 

At least within the truncation we find that the ground state occurs for = 0, 
contradicting earlier speculations and the estimates of too simple F"^ or F^ trunca- 
tions! As an immediate consequence, the gauge coupling never diverges and 

as = g'^{knp)/{^'KWi) 



-12C(3) ( N^"^ - {m^ + 8)^ + {12N^ + 23)^ ) \ (3.18) 
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reaches the value 



a,{k 0) = 3.35 



(3.20) 



4 Gluon condensate 

Let us now turn to the gluon condensate as described by the expectation value of a 
suitably smeared composite operator ~< F^^F^^ >. This expectation value may not 
vanish despite a vanishing ground state value for the gauge field < >= 0. The 
formalism for the introduction of composite fields in the flow equations is described 
in [0, and we use it here in the version without a separate infrared cutoff for the 
composite field x- Essentially, this amounts to the introduction of an identity of the 
type 

1 = const I Vxexp - J d'x^ {/cjx - \F^.Fr] (4.1) 

into the functional integral defining the effective average action r^^. The scale 
and the parameter a should be chosen conveniently such that a maximum of the 
effects of higher invariants can be described by the dynamics of the new scalar 
singlet X- Typically should be a scale where the influence of W2 and on the 
running of Wi (or g"^) becomes important. The expectation value < x > for k ^ 
is directly related |]13| to the expectation value of the composite operator F^yF^^ 

<X>= \k-' < F^^Fr > . (4.2) 

The operator appearing on the r.h.s. is regulated by smearing over distances ~ k~^. 
For scales k < k^ we now have to deal with the coupled system of and x- 
The exact flow equation has the same general structure as before, except for the 

(2) 

extension of to the second functional derivative including also the new field x- 
Nevertheless, since in the formulation used here acts only on gauge fields and 
ghosts, only the projection of (F^^-* + Rk)^^ on this restricted space enters the flow 
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equations. Formally, the components R-^^ and R^a vanish in the generalized trace 
on the r.h.s. of the flow equation. 

Including the composite field, we extend our truncation for the average action 
( |2.1|) by the addition of a term 

r,[A, x] = j d'x{V,{^F^,Fr, X) + \z^,k{xWxd,x} (4.3) 



At the scale fc^, eq. ( [4.1|) implies the "initial values" |]T3 



1 1 

_ F^^ A ( F^ F^^" 



2 



= (4.4) 

The expectation value < x > for < >= corresponds to the minimum of the 
effective scalar potential 

Uk{x) = Vk{0,x) (4.5) 

for A; ^ 0. At the scale this potential is simply quadratic, f/^^ = ^^^X^? and 
has its minimum for x = 0- We will be interested in the change of shape of Uk as k 
flows towards zero. 

As in the preceding sections we will be concerned with configurations of space- 
independent static magnetic fields where 9 = \FI^^F^^'^ is a constant. We note 
that Vk^{9,x = 0) differs from Wk^{9) by the subtraction of terms ~ 9"^, 9^ and 
9^. For all k, the original function Wk{9) can be recovered by solving the field 
equation for x, f^(^? Xo(^)) = 0, and inserting Xo(^) into V, i.e. W{9) = V{9, Xoi^))- 
Nevertheless, the introduction of x effectively extends the truncation and therefore 
results in modified flow equations for Wk{9). In particular, we can now expand 
Vk{9, x) in powers of 9 for arbitrary constant x- 
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Let us next derive the evolution equation for the scale dependence of Vk- This 
is done by evaluating the general flow equation for constant x ^md B, inserting 
on the r.h.s. the truncation ( [4 .31) . We need the inverse propagator T^^^ for this 



"background" configuration. The components (r^^'')^^^^^ remain the same as in the 
previous sections, if we simply replace Wk{0) by Vk{9,x)- The pure scalar piece is 
also easily obtained]^ 

(rf)« = i^^^ + z«.M(-a^). (4.6) 



We note, however, that the scalar fluctuations can influence the flow equation only 

k )xA^- 



(2) 

indirectly through the off diagonal piece (F^ )x-4u- Using 



we find for the constant magnetic field background 

, , f , , = {D,F^nz{x)-^^{x) = 0. (4.8) 

For this configuration we therefore find a particularly simple flow equation: Since 
(f[^^-* + Rk) is block diagonal and R^^ = 0, it reduces to the same flow equation as 
before, with Wk{9) replaced by Vfc(6',x)! 
Introducing according to ( p.l9| ) 

Vk{^,x)=9'k-''Vk{9,x) (4.9) 

the evolution equation for Vk is given by (|2.20|) , with Wk{^) replaced by ffc(^, x) 
dots denoting partial derivatives with respect to ^ at fixed x- Typically we will 
choose = knp such that g'^{k) = g^iknp). Expanding as before in cubic order in 

Vk{^, X) = u{x) + Mx)^ + lv2{x)^^ + ]:V3{X)^^ (4-10) 

2 



^°Z^(x) is related to Z^{x) and its derivative in a simple way. 
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the flow equations for vi{x),V2{x) and Vs{x) are given by ( |3.7| ) and with the 

replacement Wi Vi{x)- The evolution of u{x) follows from ( p.20| ) for i) = and 
d = 4, i.e. 

1^9\KpI_ (4.11) 



d 



dt 



u 



_4« + ^(iV2 



The system of equations (|4.11|) , ( p.7|) , and (|3.18|) still remains a complicated system 
of four coupled partial differential equations for the functions u{x,t), fi(x, t), 
V2{Xyt) and Vs{xjt) which depend on two variables. 

In order to gain some intuition about this system let us first look at the differ- 
ential equation for the scalar potential Uk{x) and its derivatives 



d_ 
di 
d_ fdU_ 

dt \ dx 
d_ f d^U 
dt \dx^ 



U 



3C(3) 



47r2 
3C(3) 



dx 



47r2 



' dvi 
dx. 



Vl 



dx^ 



(4.12) 
(4.13) 

(4.14) 



At the scale the initial values for vi and dvi / dx are given by 



dvi 



np J 



g'^ak^ 



X 



9\k. 



np ) 



g^ak^ 



(4.15) 



whereas d'^vi/dxfk^ vanishes. The r.h.s. of eq. ( [4. 13] ) starts therefore positive and 
induces a negative linear term in f/ at x = 0, resulting in a minimum at Xo(^) > 
for k < k^. Simultaneously, the positive r.h.s. of eq. ([4.14|) leads to a decrease 
of the mass term for x- One can follow the fc-dependence of the minimum value 
Xoik) ixo = Xo(^ = 0)) using the identity 



dxo 



— m. 



2 9 (dUk\ , 

(Xo 



dt ^ dt \ dx 

where we have defined the scalar mass term 

d^U 



dx- 



-iXo). 



(4.16) 



(4.17) 
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Inferring the r.h.s. from eqs. ( [4. 131 ), ( [4.14| ) evaluated at x = Xo "we conclude that 
the running of Xo stops for small k due to the factor k^, provided ml and Vi remain 
nonzero for A; ^ 0. 

Let us next consider an expansion in powers of for x at the minimum Xo(^) 
and compare 



Qn 



t;(^9,Xo(^ = 0))|^=o 



= 0) 



1^ 



v{d,xm)\i)=o 



(4.18) 



Whereas wi = vi holds for all k the corresponding relation for W2 reads 



W2 = V2 



dx^ 



iXo. 



d'd w=o' 



(4.19) 



For positive xo and negative we observe vi = fi(xo) < ^i(O). Due to the 

running Xo the flow equation for Vi obtains an additional contribution 



d g\k^^) (UN V2\ , dvi{x) d 

Vi = — — I — z -H. 



dt 



+ 



1X0 



dt 



Xo, 



(4.20) 



^ should be added on the 



and similar for all Vn = Vn{xo) where a term ^ 
r.h.s. of the flow equations. 

It is instructive to compare the flow equation for vi to the one without a con- 
densate (|3.71 ). In the vicinity of /c^ we can use the initial values (|4.17| ), together 
with 



g\k 



V2\k = ^2|fe + _4 



np I 



(4.21) 



As compared to eq. (p.7|) the difference between V2 and W2 subtracts from the 
r.h.s. of ( |4.20|) a piece g*{knp)Hi/ {^T['^g^av\). On the other hand, the contribution 
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dv 



^dxo/dt = g^{knp){H4^ — rf) / (Sn^g^avf) almost cancels this piece except for the 
negative contribution ~ rf. The flow of Wi is therefore somewhat slower than the 
one corresponding to (4.9). 

We also need flow equations for the quantities 

which appear on the r.h.s. of the flow equations for Vn, Xo and m^. They obtain by 



partial differentiation of (4.9) and (5.8) with respect to x'- 

^2/2 = 4y2 + ^^|-r2ro't;iyi + (5ri + -(iV -l)r2)(2-^-4-p) 

-{rt + liN'-lH){^-3^-^)} (4.24) 



d . g\knp)r5A8 4,3 3 174 4,2. ^ n 

+6r^(17^- + 17^^ 68 ? 105 — ^ + 175— ^ + 3^-) 

Vi Vi Vi vf 

+6r2^(iV^-l)(9^ + 9^-36^ 

-36^ + 60— g 3+3^)[ (4-25) 

As before, we truncate by neglecting terms ~ V4 and ~ 1/4 and we further omit 



d'^Vi/dx^iXo) and a term ~ ^t(Xo)^ in the flow equation for ( |4.14| ): 



1-4 = ^CV' - 1)**".''!/? (4.26) 

We then end up with a closed system of differential equations for the eight func- 
tions f 1, f2, V3, X05 2/3- This system can be solved numerically. The initial 
conditions for yi at the scale follow from (|4.4|): 

^'■».=-^- S-V.^". J/^.-" (4-27) 



g^ak^ 
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We note that the choice of is arbitrary (it concerns only the overall normalization 
of 9) and it seems convenient to take = g'^{knp)- 

In order to make some quantitative comparison we compute in addition the 
trace of the energy momentum tensor T^. In QCD its nonzero value is related to 
the running of the coupling constant and it is usually quoted as = —0.014 GeV^. 
In our context the part of induced by the gluon condensate can be expressed in 
terms of the effective potential for x, 

TH = 4[^(xo) - Um (4.28) 



Its value at = can be computed from the flow equation ( |4.12| ) 



with T^{k-x) = 0. We use a linear approximation for vi{x): i-e- ^i(x = 0) = vi—yiXo- 
As an alternative, we can estimate in the approximation of a quadratic potential 
U = |m2(x-Xo)^ i.e. 

T; = -2mlxl (4-30) 
We have solved the system of differential equations for vi, V2, v^, yi, y2, ys, Xo and 
numerically. More precisely, we have first solved for k > the system of flow 
equations of sect. 3 without composite fields. Here we have started at = 4 GeV 
with a gauge coupling adjusted in such a waythat the scale knp coincides with a 



similarly defined scale in terms of the coupling in the MS'-scheme with three light 
quark flavours: 



In practice we set 



f^a^ldiJs^knp) ^ -1.5. 



knp = 2.11 A^^j^ = 600 MeV 



(for Aqqy) = 285 MeV the two- loop confinement scale). At the scale knp we also have 
introduced the composite field according to eq. ( [4.4|) , with /c^ = knp- The results 
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depend somewhat on the parameter a that enters the initial values ( ^.151 ), ( [4. 211 ), 
([4.27|) . We consider explicitly two values of a, 

« = - -2 \ V ^ = 2(1.1) (4.31) 

where the value in brackets corresponds to a situation for which the system becomes 
numerically unstable. For A = 3 we obtain similar values as for A = 2. We consider 
the values for A = 2 as our best estimate and take the values for A = 1.1 as an 
indication of the error. 

One obtains for k ^ the same fixed point values for Vi,V2, and as as found 
before without composite operators (end of the preceding section). This does not 
change for larger values of a. For k = the expectation value xo reaches the value 

Xo = 0.65 (1.45) GeV 

for which 

<X> /Aqcd = 2.3 (5.1) 
and the "F^" condensate equals 

< F;,Fr 4fcJ <x>= [3.1 (3.7)Aqcd]' 

Here the first and the second figure always refers to A = 2 and A = 1.1, respectively. 
From the expressions ( [4.28| ), ( [4.30| ) for one finds 

Tj^ = -0.015 (-0.039) GeV^ and Tj^ = -0.0145 (-0.0295) GeV^ 

respectively. Comparing with the value = —0.014 GeV^ from QCD sum rules 
the agreement is almost perfect. Since we have made no attempt to define the 
renormalization scheme precisely and quarks are neglected, there is, however, a 
considerable scale ambiguity on top of the truncation errors. The mass term turns 
out as = 130 (83) MeV. This should, however, not yet be associated with the 
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mass of a scalar glueballQ. The latter needs knowledge about the kinetic term of the 
corresponding operator and could be defined as nigi, = m^Z~^/^. It seems possible 
to compute from the solution of its flow equation, starting at with = 0. 

We do not think that the precise numbers should be taken too seriously. Even 
within our truncation one could solve the partial differential equation for Vk{6,x) 
without a polynomial truncation. One may also investigate the dependence of the 
results on the parameters and a. This would also give a rough estimate of the 
truncation uncertainties since without truncations the final values of the condensates 
should be independent of k^ and unique for a given definition of the composite 
operator. Nevertheless, the fact that expectation values come out in a reasonable 
order of magnitude seems to indicate that our approach offers a possible perspective 
for a quantitative understanding of the gluon condensate. 

5 Conclusions and outlook 

In this paper we have approximated the exact nonperturbative evolution equation 

for Yang-Mills theories by a truncation where the effective action is given as 

a function Wk of 6* = ^F^^F^^ . More precisely, we have investigated the effective 

action for constant magnetic fields B in dependence on an infrared cutoff k and 

solved the corresponding flow equation for A; — > 0. Our ansatz is general enough 

to allow for a ground state value of B different from zero. This would correspond 

to an absolute minimum of Wk{0) for > in the limit where the infrared cutoff 

k vanishes. Polynomial approximations of Wk in order 9 and 9"^ seem to indicate 

that the minimum of Wk^^ occurs for 9q > 0, but are obviously insufficient since 

they lead to a function W{9) which is not bounded from below. The minimal 
^^The quantity should also not be confused with the correlation length for field strength 
fluctuations p^ . 
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realistic truncation in order 9^ leads to the conclusion that the minimum of Wk{0) 
stabihzes at = for all k. Prom an inspection of the general structure of the 
flow equation for Wk{9) it seems likely that this feature is not an artifact of the 
polynomial approximation but rather a property of the partial differential equation 

for Wk{e). 

For a description of gluon condensation we have, in addition, introduced a gauge 
singlet field % which is associated to a suitably regularized composite operator 
FnvF^^ . Our computation of the effective potential for x indicates a nonvanishing 
expectation value of x, related to a nonvanishing expectation value < F^iiF^'' >^ 0. 
A first numerical investigation gives a value for the condensate < F^^F^^^ >^ 
[3.1Aqcd]^ where Aqcd is the two-loop confinement scale in the MS scheme with 
three hght flavors. The precise meaning of the operator F^^F^^'^ needs to be worked 
out - in the present version it roughly corresponds to F^y6{k'^ + D'^[A\)F^^ with 

the covariant Laplacian and k^^ ~ 2.1 Aqcd- Here 9 stands for a rather sharp 
cutoff (similar to the usual ^-function) such that high eigenvalues of —D"^ do not 
contribute to the regularized operator. Even though we do not claim quantitative 
precision of our estimate we find it remarkable that the solution of relatively simple 
flow equations leads to a reasonable order of magnitude of the condensate! 

For higher precision, and, in particular, for a possible comparison with more 
phenomenological approaches such as QCD sum rules or the stochastic QCD vacuum 
several shortcomings of the present truncation should be overcome. We list a few 
which seem to us particularly important: 

(i) In the present truncation the momentum dependence of the inverse gluon 
propagator is always approximated by Z^cf, with Zj^ — WjJ(O). For k in the vicinity 
of Aqcd this is certainly not a very accurate approximation, even though the flow 
equations involve effectively only a small momentum range ~ k^ for a given scale 
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k. A more reasonable approximation for the term quadratic in the gauge field 
would be of the sort 

\jd'xF^^Z,{-D'[A])Fr (5.1) 

with Zfc depending on the covariant Laplacian in the adjoint representation^. For 
the functional form of Z].{q^) we expect for large a A;- independent positive func- 
tion Z{q^). In lowest order the logarithmic dependence of Z on should be deter- 
mined by the one-loop /3-function for the running gauge coupling. Indeed, external 
momenta of the gluons act as an independent infrared cutoff. For k"^ <^ the 
running of Zk{q^) with k should essentially stop whereas it is given by the one- loop 
/3-function for k'^^q^. In contrast, our truncation identifies Z^iq^) with a constant[^ 
Zk{q^ = 0) = Wl.{6 = 0). In the truncations of sects. 3 and 4, W^iO) remains posi- 
tive for all such that the use of Z^ = W^fc(O) induces quantitative inaccuracy, but 
no qualitative problems. 

ii) We observe that the spectrum of small fiuctuations around the constant mag- 
netic field configuration ( ^.31 ) lacks Euclidean SO (A) rotation symmetry. This is 
not surprising since F^^ singles out two of the space directions. In addition, the 
spectrum is partially discrete - continuity exists only with respect to the transversal 
momentum. The lack of full rotation symmetry and the partial discreteness of the 
spectrum are actually not unrelated: For an 5*0 (4) symmetric spectrum the conti- 
nuity in two momentum directions must extend to all momentum directions. One 
^^The exact propagator will, in addition, also receive corrections from the modified gauge-fixing 

term T gauge discussed in Appendix A. 

^•^This would lead to unphysical problems for vanishing or negative W^. With ( ^.l| ) the quantity 

Zk{Q) could turn negative without affecting the high momentum behavior of Zk{q^). With such 

a truncation the high momentum modes are always stable if the momentum dependence of the 

propagator is properly taken into account. This implies positive u for — > oo and no ultraviolet 



problem can appear for the exponential cutoff (1.2). 
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expects then a spectrum with a few separate particles, each of them having 
depending continuously on a generalized SO {4) invariant of the type g^g'^. The lack 
of full rotation symmetry for the fluctuation spectrum in a given background is not 
a problem. Nevertheless, it would be interesting to check the validity of approxi- 
mations by investigating background configurations with a different symmetry. In 
ref. [0 we have proposed a candidate for a configuration with generalized 5*0(4) 
symmetry for d = 4 and gauge group SU{N), N > A. Unfortunately, no realistic 
candidate for the gauge group SU (3) has been found up to now. It would be very 
interesting to perform an analysis similar to the one presented here for SU{A), using 
the SO {4) symmetric configuration instead of ( p.3|) . 

For the gauge group SU (3) an interesting alternative configuration with constant 
colour-magnetic field is given by 



(with i = 1, 2, 3 corresponding to the spacelike indices in Minkowski space). It is in- 
variant under the diagonal subgroup 5*0(3) of three-dimensional space rotations and 
global gauge transformations (with generators Ti, T2, T3 forming an SU{2) subgroup 
of SU {?>)). The field strength corresponds to a constant colour- magnetic field 



adf for z = 1, 2, 3 



(5.2) 



for 



I > 3 




for i,j = 1,2,3 



= < 







otherwise 
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Fr^;, = 2B' = Q-g'a^ (5.3) 
Nevertheless, F^y is not covariantly constant {i,j,k = 1,2,3) 

iF''.,,). = -fa'm-6i6l) (5.4) 

i-D'[A]rF:^ = ^fa'py^ (5.5) 

and this configuration allows to explore a nontrivial momentum dependence of the 
gluon propagator as in ( ^.1| ). For the characteristic properties of the spectrum of 
fiuctuations around the configuration ( |5.2| ) it is sufficient to restrict the discussion 
to = 3. With the truncation ( p.l| ) we find that the spectrum is continuous. For 



large momenta it turns out positive semi-definite for all B, even for k = 0. This 
contrasts with the spectrum for the background ( p.3|) ! For a three-dimensional the- 
ory (high temperature field theory) the configuration ( ^.2|) is an interesting ground 
state candidate. For the four- dimensional theory a study of Tk as a function of a 
should reveal additional aspects of the functional form of Ffc. 

In general, an investigation of the functional dependence of F^ on various distinct 
configurations would greatly enhance the robustness of our results. Qualitative and 
quantitative results should be independent of the configuration chosen. The use 
of different configurations should therefore permit an estimate of the truncation 
uncertainties. In addition, different configurations project on different invariants 
and the use of several configurations would give much more detailed information 
about the full functional form of the effective action F^. 

iii) The present treatment of composite operators is still rather rough. One would 
prefer to work directly with a smooth definition of a regularized operator, as, for 
example, 

X~F,yexpiD'[A]/kl)F^'^ (5.6) 

X 
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This would permit a more direct comparison of the expectation value of x with a 
corresponding weighted momentum integral of the two-point function for F^^, and 
therefore with lattice results or phenomenological estimates. The general formalism 
for the treatment of this type of operators is known [|l^], but in practice it will 
require a more detailed study of the momentum dependence of F^, beyond the 
truncation ~ Wk{0). In fact, for the configuration with constant B used in this 
paper D'^[A\F^'^ vanishes. At the present stage we can therefore not distinguish 
between different versions of a regularized F^j^F^^'^ operator. (This is the reason why 
we cannot specify the function 9 mentioned above.) In order to resolve this issue, 
one needs an investigation of different configurations as described in ii). 

In summary, our first attempt to investigate the phenomenon of gluon conden- 
sation with the help of nonperturbative flow equations is encouraging. The simple 
configuration ( p.3|) and the truncation (2J.) give a qualitatively interesting picture: 



Whereas the ground-state value of the gauge field vanishes, a composite operator 
~ F^iyF'^'^ develops a non- vanishing vacuum expectation value. For the moment our 
analysis is too rough for a detailed comparison with quantities that can be deter- 
mined in lattice QCD simulations or by phenomenological estimates as QCD sum 
rules. We hope that a future extension of our investigation will allow for such a com- 
parison. We find it remarkable that already a relatively simple scheme which starts 
from first principles and leads to a quantitative estimate of the gluon condensate. 

Acknowledgement: We would like to thank D. Jungnickel for a first numerical 
solution of the differential equations of sect. 3. 
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Appendix A 

In this appendix we derive the exact form of the evolution equation for as well 
as the modified Ward-Takahashi or Slavnov- Taylor identities and "background field 
identities" which it satisfies. We start from the scale dependent generating functional 

exp Wk[K;, a\ a'- 7^; A^] = J VAVCVC exp - {S[A] + A^S 

} = Jvc^exp-Stot (A.l) 

where S[A] denotes the gauge invariant classical action, and 

A.^ = lJd''xiA-A)lR,{A)ll{A-A): 

+ J d'^xcyRk{A)y'C' (A.2) 

is the infrared cutoff for the gauge field fluctuation a = A — A and for the Faddeev- 
Popov ghosts C and C Here Rk{A) is a suitable cutoff operator which depends on A 
only. It may be chosen differently for the gauge field and for the ghosts. Furthermore 

^gf = ^ / d'x [D,iA)y^{A - A);] ' (A.3) 

is the background gauge fixing term and 

^ghct = - / d'^xCy {D^{A)D,{A)y' (A.4) 

is the corresponding ghost action [|l^]. The fields A — A, C and C are coupled to 
the sources K, a and a, respectively: 

5'source = - \ d" x{KI{A'^ - A^) ^ ^ 

+^^yD^{AY'C' + ^^^f^y^cyc'}. (a.s) 

(In this appendix we use throughout g instead of g.) We also included the sources 
P and 7 which couple to the BRS- variations of A and of C, respectively. In fact, 
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^ + ^gf + 'S'ghost is invariant under the BRS transformation 

2 ^ 



6C- = ^D,{Ar{A;-A;). (A.6) 
Let us introduce the classical fields 

and let us formally solve the relations a = a{K, a, a; P,j;A), ^ = ^(...), etc., for the 
sources K, a and a : K — K{a, ^, ^; ^,^;A), a — .... We introduce the new 

functional Ffe as the Legendre transform of Wk with respect to a and a: 

f,[a,e,e;^,7;^] = Jd'^xiK^a^ + aT + ^T] 

-Wk[K,a,a-J,r,A]. (A.8) 

Apart from the usual relations 
we have also 



5tk 5Wk 5Vu 5Wk 



(A.IO) 



where 5T/5(3 is taken for fixed a, ^, ^ and for fixed cr, etc. 

The effective average action is obtained by subtracting the IR cutoff /^kS, 
expressed in terms of the classical fields, from the Legendre transform r^: 

-J d''x^yRk{A)y'C- (A.ii) 
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Frequently we shall use the field A = A + a (the classical counterpart of A = A + a) 
and write correspondingly 



Tk[A, A, e, e; P, 7] = - A, e, e; p, % a]. 



(A.12) 



For ^ = ^ = /? = 7 = one recovers the effective average action Tk[A, A] introduced 
in ref. [0. Using the methods described there one finds the following exact evolution 
equation governing its scale-dependence, with t = \nk 



^Tk[AA,^,^-J,^] = ^TT,,L 



{r^^^ + R,{A)y^\^^R,{A)AA 



1 

Tr 

2 



(A. 13) 



(2) ^ ^ ^ — 

Here is the Hessian of F^ with respect to A, ^ and ^ at fixed A, j3 and 7 and 
RkAA,Rk^^ are the infrared cutoffs introduced in ( |A.2|) . (We use for the ghosts the 



convention (f|^)^^ = ^f^, (f[|) 



•) 



It is clear from its construction that F^ is invariant under simultaneous gauge 
transformations of A^ and and homogeneous transformations of ^,^,(3^ and 7, 
i.e., 5Ffc[A,Ae,e;A7]=0for 



6vy 



--^D^{A)y'uj' 
-^D^{A)y'uj^ 



(A.14) 



Next we turn to the Ward identities. By applying the transformations ( |A.6|) to 
the integrand of ( |A.1|) one obtains from the BRS invariance of the measure Vcj) 



I V(I)5brs exp -Stot = (A. 15) 
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or 



/ d'x Ik; 



sw. 



(fx 



.SWk 



+ 



+ 



H 

6a^ 6a 



5 



R, 



^7 



SWk s 

+ 



6K„ 6K„ 



Rk 



6Wu 



6a J 
(A.16) 



with {Rk6Wk/6^y = Rk{Ayy6Wk/6^y, etc. Equation (^]Tg) can be converted to 
the following relation for the effective average action ( |A.12| ): 



A 



(BRS) 



where the symmetry-breaking contribution A^,^^^"* is given by 



(A.17) 



A 



(BRS) 



-Tr 

-^Tr 
ag 



2-p/ 



-1 S'T 



RkiA),,{Tf^ + Rk,^^^^^^ 
D,{A) {Tf+Rk)~^ Rk{A)-^^ 



and where 



V'k^Vk--Jd^x 



D,{A){A,-A,] 



(A. 18) 



(A.19) 



Here ip = (A^,^,^) is summed over on the r.h.s. of ( [A.18| ) and the traces on the 
r.h.s. of (A. 18) involve an x integration as well as a sum over the suppressed index 
z. In deriving eq. (A. 18) we used 



Tk[A,A,^,^;P,^] = 



(A.20) 



which follows from the equation of motion of the antighost. 

Equation ( [A.IT] ) is the generating relation for the modified Ward identities which 
we wanted to derive. In conventional Yang-Mills theory, without IR-cutoff, the r.h.s. 
of (|A.17|) is zero. The traces on the r.h.s. of (|A.17|) lead to a violation of the usual 



38 



Ward identities for nonvanishing values of k. As k approaches zero, Rk and hence 
A^,^^^'' vanishes and we recover the conventional Ward-Takahashi identities. Eq. 
(|A.17|) is equivalent to a similar identity derived in ref. |ll|] using a different gauge. 



The modified Ward identities (|A.17| ) are not the only conditions which the av- 
erage action has to satisfy. Next we derive a relation which constrains its de- 
pendence on the background gauge field A^. Taking the A-derivative of (|A.1|) one 
obtains 



-{Rk-aYJy) + -{D ® Dayjy) + Kiy) 



a 



1 6{Rk-^D^D) 



SAM 



-Or +^ 



SR, 



^ (Rk --D®D) 



a 



-Tr, 



5^Wk 6 



5a 5 a 5Ai{y) 



Rk 



-9Tiy)f 



vzw 2 fzvw 



(A.21) 



6(3^{y) 6a\y)6(3^{y) 
with {D ® Dyp = D^iAyW^iAy". Upon Legendre transforming eq. and 
switching from to F^, one arrives at the following result 



-Tr 



6 



AA SA^^iy) 

1 



I (fx tr^ 



t^{t^^^+rx: 



Rk--D®D 

OL / AA 



(A.22) 



Again, (p = (A,^,^) is summed over and (T^)^^ = —if^y^. Note that the r.h.s. of 
eq. ( |A.22| ) does not vanish even for ^ 0. The D ® Z)-piece of the 2nd term and 
the 4th term on the r.h.s. of (|A.22|) survive this limit. 
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So far we were deriving general identities which constrain the form of the exact 
functional F^. Let us now ask what they imply if we truncate the space of actions. 
In the present paper we neglect the fc-evolution of the ghost sector by making an 
ansatz which keeps the classical form of the corresponding terms in the action: 

TkiA, A, e, e; 7] = ^k[A, A] + Tgh (A.23) 

Tgh = - / d'^x ^D^{A)D^{A)^ 

-J d'^xi^^^J^{Df,{Ay'"C + 17"/™'"^^^"'! (A.24) 

If we insert this truncation into the general evolution equation (|A.13|) , we obtain 
precisely eq. ( |l.lj ), which is our starting point in the main body of the paper. A 
generic functional rfc[y4,yl] can be decomposed according to 

T,[A, A] = T,[A] + ^J d'x[D,{A)iA, - A,)]' + f™[A, A] (A.25) 

where is defined by equating the two-gauge fields: Tk[A] = rk[A,A]. The re- 
mainder r/c[y4, y4] — Tk[A] is further decomposed in the classical gauge-fixing term 
plus a correction to it, ff'^"^'^, which also contains the counterterms. Note that 
-pgaugcj-^^ _ g equal gauge and background fields. We observe that Tk[A] is 
a gauge-invariant functional of A^ and rfc[yl,A] is invariant under a simultaneous 
gauge transformation of A and A. 

In the present paper we make the further approximation of neglecting quan- 
tum corrections to the gauge fixing term and fc-dependent counterterms by setting 
pgauge ^ Q_ rj.^^^ ( |I35| ) coincides with eq. (^) for 



Tk[A] = I d^x Wki^F''). (A.26) 

The important question is whether this truncation is consistent with the Ward- 
Takahashi identities ([A.17|) and the A-derivative ( [A.22|) , respectively. If we insert 
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( [A.23D -( |A.25[ ) into ( [A.!?] ), for instance, we find that drops out from the l.h.s. of 



this equation. We are left with a condition for f l'^'^^'^ 



1 /■^d^^rs^^s'^ 



(BRS) 



(A.27) 



The term A^,^^^"* ([A.18|) vanishes for — » but is non-zero for A; > 0. Our approx- 
imation f l*^"^*^ = is consistent provided these terms can be neglected. We note 
that the traces appearing in ( |A.27| ) are related to higher loop effects. Beyond a loop 
approximation our neglection of r^'^"^'^ is a non-trivial assumption. We emphasize 
that because of its gauge invariance the functional Tk[A] does not appear on the 
l.h.s. of the Ward identities. Therefore the Ward identities do not imply any further 
condition for r^. This means that, within the approximations made, we may write 
down any ansatz for Tk as long as it is gauge-invariant. 

Similar remarks apply to the identity for the A-dependence ( |A.22| ). If we in- 
sert ( |A.23| ) with (|A.25|) into (|A.22| ) the first term on the r.h.s. is cancelled by 
6Tgi^/6Ajj^{y). One obtains 



gauge 



SAliy) 



2a 



+g^ J d'^x trc 
+Tr[Ru{...)] 



^{yMx) 6(p{x)6/3^{y) 



(A.28) 



Consistency of the truncation r|^"^^ = requires that we neglect the traces on 
the r.h.s. of ( |A.28|) . Contrary to the case of the Ward identities, not all of these 
terms vanish obviously for k 0. We observe again that, within the present 
approximation, the equation for the A-dependence does not impose any restriction 
on f fc. 

We finally observe that the flow equation ([lT|) can be rewritten in close analogy 
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to a one-loop formula: 



d_ 

dt 



2Dt 
D 



1 D 



■Tr,,ln \-D''[A]D^[A] + Rk{-D^[A]) 



Dt 



(A.29) 



The derivative jyj: acts only on the explicit ^-dependence of the function Rk, but not 
on r^,^"* [A, A] . It is now easy to describe the relation between the effective average 
action and the conventional perturbative effective action. Let us first briefly 
discuss the approximation — > ^ in eq. (|1.3| ). This amounts to neglecting the 
running of on the r.h.s. of the evolution equation. It is then trivial to solve it 
explicitly: 



with A some appropriate high momentum scale (ultraviolet cutoff) where we may 
identify Fa with the classical action S including a gauge fixing term and coun- 
terterms. This formula has a similar structure gularized expression for the 



two important differences, however: 

(i) The second variation of the classical action, 5"*^^^ is replaced by F^^'*. This 
implements a kind of "renormalization group improvement" and transforms ( |A.29|) 
into a sort of "gap equation" . 

(ii) The effective average action contains an explicit infrared cutoff Rk- For the 
choice ( |1.2|) one has 




conventional one- loop effective action in the background gauge [^, 0]. There are 



lim Rk{u) = 0, 



hm Rk{u) = Zkk'^. 



(A.31) 
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(2) 

Effectively, a A;-dependent mass-type term is added to the inverse propagator 
for tlie low frequency modes (m — > 0), but it is absent for the high frequency modes 
{u —>■ oo). Despite the similarity of (|1.3| ) with a one-loop expression, we stress that, 
for /c — > 0, the solution of the original renormalization group equation where D/Dt 

(2) 

does not act on equals the exact effective action which includes contributions 
from all orders of the loop expansion^^. 

For a detailed discussion of the approximation ( |A.29D in the case of the abelian 



Higgs model we refer to , and to ref . [|l5l for the corresponding nonperturbative 
evolution equations of this model. 



Appendix B 

In this appendix we discuss the technical steps needed for the derivation of the 
truncated flow equation ( |2.5| ) from the exact equation ( |A.13| ). Upon performing the 
second variation of the ansatz (p. 11), 



5'T,[A, A] = l d'^x5A'^T^^\A, A^^dA^ (B.l) 

(2) 

we arrive at the following operator F^ : 



Ff [A, A]ll = W;,{9){Vt[A] - V,[A])y;^ + WmSll[A] + Uvl[A])11 (B.2) 



with 



= \F^.Fr (B.3) 



Here we used the notation {w, y, z are adjoint group indices and g is the (bare) gauge 
coupling) 



^''This holds for the exact flow equation (|A.13 ), whereas (1.1) involves already an approximation 
in the ghost sector. 
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{Vl)% = -iD0D)ll = -Dy-Dr 

Sf. = Fy^FZiDWT' (B.4) 

with the covariant derivative {Dn[A]y^ = 8^6^^ — igA^^iTzY^ in the adjoint rep- 
resentation and = F^j^iTzy^ . Moreover, and W^l denote the first and the 
second derivative of Wk with respect to 6. 

In writing down eq. ( |B.4| ) we made the additional assumption that the field 
strength is covariantly constant F^^.^ = or 

[Dp[AlF^,[A]]=() (B.5) 



It is easy to see that ( p. 3D and (|2.4| ) obey the condition (|B.5|) . The choice (p.3| ) has 
the advantage that it allows for an explicit diagonalization of the operator . 

We note that can be extracted from rfc[yl,y4] which is a gauge-invariant 
functional of A obtained by putting A = A. It is therefore sufficient to know 

Before turning to the evolution equation, we list a few special properties of the 
covariantly constant fields, which will prove helpful later on. From (|2.3|) it follows 
that Afj^ satisfies the classical Yang-Mills equations D^F^y = 0. This in turn is 
sufficient to prove that the operators Vt and commute. As a consequence, one 
may define generalized projection operators [0 

Pl = V^'Vl 

Pt = 1-Pl (B.6) 

which satisfy P^^ = Pt,l, Pt + Pl = I and PtPl = = PlPt- For = they 
reduce to the standard projectors on transverse and longitudinal modes: 

{pP)^.u = s^,-d,d,/d^ 

(pPU = dMd' (B.7) 
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Furthermore, if is of the form ( ^.3| ), it is natural to define another pair of orthog- 
onal projectors, 

pyz ^ ^yz _ ^y^z^ pyz ^ ^y^z^ ^g g^, 

which project on the spaces perpendicular and parallel to n^, respectively. For the 
vector potential (2^) the matrix A^T^ reads in the adjoint representation 

AJ^ix) = {A;Tjy^ = z/^"%^A^(x) (B.9) 

The antisymmetry of the structure constants f^^^ implies that P\\ and P± commute 
with D^, D^, Vt, T>l and F^^,, and that 

P||A^ = 0, P\\D^ = P\\d^, P\\{VT)^u = -d^^^.P\\ (B.IO) 

The operator S from ( p.4|) factorizes according to 

cyz _ pyz 
'-'^u — -'II ^fiiyj 

s,u = f.pfaud'd'^ (B.ll) 

Hence S commutes with D^,Vl and Vt because it annihilates the gauge-field con- 
tained in these operators: 

SD^ = Sd'^, SVT = -Sd'^, SD(^D = Sd(^d (B.12) 

In physical terms this means that those components of the gauge fluctuations 
6A^^ = which are parallel to decouple from to some extent. In fact, in 
terms of the projections a^'" = P±,\\a^ the quadratic action ( |B.1|) with ( |B.5|) reads 



5'Tk[A, A]= f d^x \af[-d^Wl5^, + {W^ - —)d^d, + W'^t 
J I ak 

+ai^mVT + (- - W;,)V,]llaj'^] (B.13) 

Oik J 

We observe that the a"-modes couple to the external field only via the derivatives of 
Wk = Wk^i^B"^). In a conventional one-loop calculation one uses the classical Yang- 
Mills Lagrangian ^-F^^, rather than Wki^F"^^). In that case the quadratic action for 
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the small fluctuations is given by ( |B.13| ) with W^. = 1 and Wj! = 0. Hence the one- 
loop determinant resulting from the integration over a" is field- independent and may 
be ignored. In the present case, the a"-modes are important for the "renormalization 
group improvement", however. 

The quadratic form ( |B.13| ) can be diagonalized even further by introducing the 
longitudinal and transversal projections 

all'^ = P^all = Pf all, a^'^ = Pra^ = P^'^^ 
a^'^ = Plq^, a^'^ = Pra^ (B.14) 
By virtue of = sPf ^ = 0, P^ s = sP^ = s and [Pl(t), ^] = one obtains 
6^rk[A,A] = Jd''x{ay-''[-d^WXu + W;!s, 



+ (B.15) 

This block-diagonal form of V), will facilitate the evolution of the traces occurring 
in the evaluation equation. For example, a"'^ gives no A-dependent contribution 
and, except for an irrelevant constant, the only dependence of on arises from 
a-^'^. Writing 

F(') [A, A] = Ff ) + Ff + + Ff (B.16) 
where, in an obvious notation 

Ff) = P||PTF[|'^i^|PT, F^) = /^|PLF[|'[P||Pi 

= PxPtFIlV^Pt, Ff = P±PlT^I\p±Pl (B.17) 
with [P\\,±,Pl,t] = and 

rf Fg) = for A^B 

[Ff,Fg)] = (B.18) 
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one obtains ([r||^^, /^iPy] = etc.) 

i?fc(r(2)) = P\\PTRk{r\;}r)P\\PT + PiiPLRk{r\^l)P\\PL 

+P^PTRk{T%)P^PT + PLPLRk{Tl\)P±PL (B.19) 



We will choose the matrix Zk in the definition of Rk ( |1.2| ) as 2^ = 1 for the 
ghosts and 

Zk = ZkPT[A] + ZkPL[A] (B.20) 

for the gauge boson degrees of freedom. Here Zk-, Zk are fc-dependent constants 
and we observe that the choice ( [B.20| ) is compatible with ( [B.19| ). If we insert the 



truncation (|2l] ) into (iLlf) with A = A,we obtain (6 = \B^) 



^^Wk{e) = ^-TT,,,[H{r^\A,A])] 

-Tt^,[Hg{-D^[A])] 



+ ^Tr.cL[PxPL(^(rit) - H{r^ll))] (B.21) 



where is given by ( [B.5| ), H{u) is defined by eq. ( PTT] ) and ^2 = / d'^x. 

Let us pause here for a moment and derive a set of trace identities which will be 
needed for the evaluation of ( |B.21j ). For the covariantly constant background ( |2.3| ), 
eq. ( [B.IUD implies for any function / 

Tr,,z.[P||/(D^,P||,Px)] =Tr,i[/(9^,l,0)] (B.22) 

because Trc[P||] = n^n^ = 1. Writing Pj_ = 1 — Py and exploiting S oc P\\ it is also 
easy to see that 

TT,cL[P±f{D^,S)] = Tr,eL[/P^,0)] - TT,L[f{d„0)]. (B.23) 

>(0) 



Since F^,y is a constant matrix, the operator of (|B.11D commutes with Pi and 



Pj> and satisfies d^s^u = 0. This fact can be used to show that 



>(0) 

,{0) p(0) 



Tr.^[/(Pr,Pr;^)] 
47 



= Tr.L[/(0, 1; s)] + Tr.[/(1, 0; 0)] - Tr.[/(0, 1; 0)] 
If one subtracts the same expression with s = one obtains 

= Tr,i[/(0,l;s)]-dTr,,[/(0,l;0)] 



(B.24) 



(B.25) 



In the last step we used that Tt^l = dTr^ for an operator ~ S^^. In the above 
identities the function / may also depend on further operators provided they com- 
mute with those displayed explicitly and do not introduce any additional colour or 
Lorentz index structures. 



For the evaluation of Ui in eq. (|B.32D we need another important relation: 

TT.cL[PLf{VT)] = Tr,,[/(-D2)] (B.26) 

It follows from the fact that the operator {VT)f_iu = —D^^nu + '^igPfiu, when restricted 
to the space of longitudinal modes (a^ = {PL)'^a^), has the same spectrum as —D^ 
acting on Lorentz scalars. The proof makes essential use of the identity 



VriD ®D) = {D® D)Vt = -{D ® D){D ® D) 



(B.27) 



which holds true whenever the gauge field contained in the covariant derivatives 
obeys D^'F^, = 0. 

Equipped with the above trace identities, we now resume the evaluation of the 
flow equation ( |B.21| ). The first trace on the r.h.s. of eq. ( [B.21| ) can be simplified as 
follows. Inserting a factor of 1 = P\\ + Pi leads to the decomposition 

'H{Tf[A,A])\=Tl+Tt 



Tr, 



xcL 



;b.28) 



with 



Tr^L 



P\\H{WIJ)t + [Wl ]D ®D + W'^S) 

H{-Wld'' + [W'k - —]d ®d + W'^s 



[B.29) 
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where ( [B.22| ) was used, and 

-Tr^.L 

where ( [B.23| ) was exploited. Let us write 



P^H{WIPt + [Wl ]D ®D + W'^S) 

H{WIVt + [Wl - —]D ® D 

H{-wid^ + [wi - —]d ® d) 



(B.30) 



Tr 



xcL 



(B.31^ 



with Ui the "nonabehan" trace 



H{WIVt + [Wl ]D ® D) 

Oik 



Ul — Tr a;cL 

and U2 the sum of t| and the second term of ( p.30|) : 



(B.32) 



Tr^L 
-Tr^;!, 



H{-d'[w;,pp + -pn + w;:s 



(0)i 



Oik 



(0)- 



:B.33) 



It is quite remarkable that if we now apply the identity (p.25| ) to t/2, the longitudinal 
contribution drops out completely and the result becomes independent of the gauge 
fixing parameter a^: 



U2 = i:T,L[H{-d'wl + w'^s)] - di:i,[H{-d'wl)] 



(B.34) 



The operators entering (|B.34|) are easily diagonalized in a plane-wave basis. A 
standard calculation yields, for W'^ > 0, IV^ + B'^WH > 0, 

-1 

/o 



1 



/ dx x^-^Hix) (B.35) 

JO 



with Vd = [2''+^7r'^/^r((i/2)]^-'^. (As always, the argument of Wk and its derivatives is 
understood to be 
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Next let us simplify the trace Ui by inserting a pair of projectors: 



Tr 



xcL 



PtH{Vt[W;,Pt + —Pl]) 



+ Tr 



xcL 



PlH{Vt[WIPt + —Pl\) 



(B.36) 



The a-dependence of Ui is contained in 

A[/i(afc) = Tr^cL 
= Tr 



1 



Pl{H{—Vt)-H{W'J)t)] 
ak 

1 

k 



H(-—D'^)-H(-WiD^ 



ak 



(B.37) 



In the last line of (|B.37|) we made use of the identity (p.26| ). 

By a similar combination of the trace identities we can also evaluate the last 
term on the r.h.s. of ( |B.21| ) 



I \ \ak J \ak 

( - H ( -^^]\ 



-Trx H 



ak 



H 



ak , 



[ \ ak/ \ «fc / j 
Evaluating the second term in a plane wave basis yields 



[B.38) 



^Tr. I H (-^ - H (-^^ 



Vd dxx^-UH( — ]-H 
Jo \ \akJ \ak 



X 



X 



[B.39) 



At this point we have exploited the various trace identities as much as possible. 
Combining these results yields the flow equation 



Appendix C 

In this appendix we discuss the group theoretical factors J2e ^j"^ appearing in the 
Euler-McLaurin expansion of the spectral sums. 
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The LHS of the evolution equation is dtWk- The argument of Wk is \F^j^F^" = 
^B"^ which is manifestly independent of the unit vector which specifies the di- 
rection of the field in "color space". The r.h.s. of the evolution equation consists 
of expansions such as which involve the factors J2e^e"^- As {ui} are the 



eigenvalues of n^T^, we can rewrite them as 

^ i.2m ^ . . . ^...Trjr^.T,, ■ ■ ■ T,,J (C.l) 

e 

where the trace is in the adjoint representation. The question is whether the in- 
variants (|C.1D are all independent of the direction of n^. In appendix D we explain 
in detail that generically this is not the case. If the orbit space of the gauge group 
in the adjoint representation is nontrivial, different n's can lead to different sums 
J^e^i"^- The resolution to this puzzle is as follows. For m = 1 we can use the 
standard orthogonality relation 

to prove that J2e = Nn^riz = N is independent of the direction of n. Likewise, if 
the symmetric invariant tensor Trc T(^zi ■ ■ ■ T,,^^) is proportional to the trivial one, 
S{ziZ2^zsz4, ■ ■ - ^Z2m-iz2,n)y ^au agalu use the normalization condition n^riz = 1 to 
show that the r.h.s. of (|C.1|) is independent of n. 

The situation changes if there exists a totally symmetric invariant tensor %iZ2---z2m 
which is different from the trivial one. Then we might have 



Tr. 



'Tm^{ziZ2 ' ' ' ^Z2m-lZ2m) '^Z\Z2- ■ ■ Z2m, (^-3) 



with some coefficient r^. (If there exists more than one T an appropriate sum is 
implied.) In general n^^n^2...'7^i22- - ^ill be direction dependent [^. If some invariant 
tensor T exists, the correct way of deriving the evolution equation for Wk is to 
compare coefficients of a fixed tensor structure on both sides of the equation. Clearly 
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the l.h.s., dtWk{\F^j^) , gives rise to the trivial tensor structure only. Therefore only 
the Tm-piece of ( |C.3D should be kept in (|C.1|) and the (n-dependent) part coming 
from T has to be discarded. Thus may be used on the r.h.s. of the equation 

for Wk provided we interpret J2 as the coefficient Tm- 

On the other hand, a nontrivial T permits us to construct additional invariants 
from an even number of covariantly conserved F^j^. Then the truncation W{d) is 
not sufficient any more to parametrize the most general effective action for constant 
magnetic fields of the type introduced in sect. 2 (with covariantly constant F^y). 
The evolution equation for the new invariants can now be extracted by projecting 
the r.h.s. on the appropriate tensor structure. We will not pursue this generalization 
in the present paper. 

In appendix D we show that for SU{2) this complication is absent. There exists 
no additional invariant tensor T, and one finds the n-independent result 

E^r = 2 (C.4) 



1=1 



for all m = 1, 2.... 



Appendix D 

In this appendix we investigate in more detail the group-theoretical quantities J2i ^'f™ 
which occur in many calculations involving covariantly constant backgrounds of the 
type = n^k^. Here we consider an arbitrary (semi-simple, compact) gauge group 
G with structure constants f^y^. For a fixed unit vector we consider the matrix 

^yz ^ n'"{T'"y' = ify'^^rf (D.l) 

The numbers z/;,/ = l,...,dim G are the eigenvalues of n: n^^tpf = uiipf . This 
equation can be rewritten in a more suggestive form. Let denote the generators 
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of G in an arbitrary representation: [t"',/!;^] = if^^^t^. If we define 



(D.2) 



the eigenvalue equation becomes 



[n, ipi] = uiipi 



(D.3) 



Clearly the z//'s do not depend on the representation chosen. We would like to know 
how the spectrum {ui} depends on the vector n. First of all, it is clear that if V 
is any group element in the ^-representation, the matrices n and n' = VhV~^ have 
the same spectrum, i.e. the spectrum is constant along the orbit of G in the adjoint 
representation. If two directions n and n' are not related by a group transformation, 
then the spectra {vi{n)} and {vi{n')} can be different. Typically, for G large enough 
[ P^ , the orbit space is indeed nontrivial, and the spectrum "feels" the direction of 



Let us go over from the basis {T^} to the Cartan-Weyl basis {iJj, i?^} of the 
abstract Lie algebra. Here a G R'^ are the root vectors and i = 1, = rank G. 
For definiteness we assume that the t^'s are in the fundamental representation where 
we write {hi, e^} for the Cartan-Weyl basis. Thus 



We assume that the generators hi of the Cartan subalgebra are given by diagonal 
matrices. By an appropriate transformation n VnV~^ any n can be brought to 
diagonal form. Therefore, in order to investigate the n-dependence of {vi{n)}, it is 
sufficient to consider fi's which are in the Cartan subalgebra: h = I]i=i ^ihi. For 
this choice 




aiE^ and [hi,e^] 



(D.4) 




(D.5) 



53 



and the nonvanishing eigenvalues ui = are given by 

r 

i^d = '^niai (D.6) 

i=l 

Therefore the quantities J2i ^f"^ can be computed exphcitly from the root system: 

dimG (rank G \ ^"^ 

E E E (D.7) 

1=1 Toots{a} \ i=l / 

Let us consider a few simple examples. For G = SU{2) we have r = l,ni = 1 
and there are only two (one-component) roots: a = ±1. Thus 

^j.f'" = 2, m= 1,2,3,... (D.8) 
1=1 

depends neither on the direction nor on the power m. This degeneracy can be 
understood by noting that for SU (2) the square of the matrix ( p.l|) is the projector 



Pr. nn = P±. This means that Ei^f"" = Tr(Pf) = Tr(P^) = 2, as it should be. 
Contracting eq. (|C.3| ) with n^^...n^'^"^ and comparing the result to ( |D.^ ) we see that 



there exists no nontrivial invariant tensor T and that = 2 for all m. 

For G = SU{3) we have r = 2, and a 2-component unit vector (ni,n2) specifies 
the direction of the field in the Cartan subalgebra. Using the explicit form of the 
roots it is straightforward to derive that 

^ zyf- = 2^-2- [{n, + 73^2)'"^ + [n, - VSn^)^^ + {2n,f^] (D.9) 
1=1 

For m = 1 and m = 2 it turns out that this expression depends on rii and n2 only 
via + n| = 1, and one obtains the direction-independent results 

E^f = 3, E^f = | (D.IO) 

1=1 1=1 ^ 

Starting from m = 3, the invariants are explicitly n-dependent. Writing ni = 
cos 9, n2 = sin 6 we find for m = 3 



8 Q r 

5^ z/f = — 11 cos^ + 15 cos^ 9 sin^ + 45 cos^ sin^ + 9 sin^ 



16 



fD.ir 
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As discussed in Appendix C, the n-dependence is related to the existence of 
a nontrivial invariant tensor T^^ za- However, we are not going to calculate the 
corresponding coefficient ts here. On the other side we note that the definition of 
T in ( |C.3| ) depends on the specific truncation. Different truncations - i.e. different 
choices of the definition of a term ~ %i...z6F^lF^lFp^Fpi^F^lF^l whose coefficient is 
set to zero - correspond to a different angle 6 for which the r.h.s. of ( |D.11| ) equals 



T3. For = 3 we may therefore use r2 = | and rs between || and If, the last 
uncertainty reflecting the uncertainty from this particular part of the truncation. 

Appendix E 



In the regime where gB/k^ ^ 1 the use of the Euler-McLaurin expansion ( |2.11| ) 



becomes questionable and we should look for an alternative representation of the 



spectral sums (|2.9| ) and ( p. 101) . In this section we use the Schwinger proper-time 



representation |^0|- It can be easily applied only for the simplifled cutoff function 

Rk{x) = Zkk\ (E.l) 
In this case one may write {x = T>t) 

Inserting this representation into (|2.9|) we may employ 



00 2 
Y,exp[-sW;,g\ui\B{2n+l)] 



n=0 

and 



2smh[sWl^g\iyi\B] 



Jo 

as long as Wj. > 0. One flnds 
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d ,„ , , f°°ds (4-d)/2 -sZkk^ 



^ exp{-sWl,g\ui\B) + exp{+sWl,g\ui\B) 



(E.3) 



2smh{sW;^g\ui\B) 

where the last exponential is due to the unstable mode. For Z^k"^ < WIj)\vi\B 
it makes the s-integration divergent at the upper (i.e. IR) limit. In conventional 
calculations of the one- loop effective action this creates a problem from the outset, 
because one attempts to work at fc^ — ^ there. In the present formulation everything 
is well defined for /c^ sufficiently large, and one interesting question is how the 
renormalization group flow behaves as one approaches k"^ 9 from above. 

In the UV limit s ^ the terms inside the square bracket in ( |E..j| ) behave as 

' ' - + 0{s) (E.4) 



2-g\ui\BWl{B^l2) s 

While the 0(s)-terms do not lead to UV-divergences for d < 6, the term ~ 1/s leads 
to a divergent contribution to the proper-time integral. Though the factor gli^ilB 
cancels against a similar one coming from the density of states, this divergent piece 
is still field-dependent because of the -B-dependence of Wj^. This UV divergence 
shows a failure of the truncation for the mass- type cutoff function Rk = Z^k"^. The 
latter may be used only together with the approximation W^l = on the r.h.s. of 
the flow equation. The divergent piece in the proper-time integral is an irrelevant 
constant then. 

Using ([E.3|) and a similar formula for the scalar traces in eq. ( p.5|) , we obtain 
(up to an irrelevant constant and for Z^ = l/a^jfj = 0) 



{2-r^)Z,{WlY-i 



^ ^ +2s\-ah{sWlg\ui\B) 



2sinh(sVr^^|z/;|S) 



exp(— sA;^ 



smh.{sg\vi\B) 



(E.5) 
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This evolution equation is the analogue of ( p^.l3D with the additional assumption 



Wj! = 0. Contrary to the Euler-McLaurin series it is valid even for strong fields 
gB ^ k"^. 

For gB ^ the r.h.s. of ( [E.5| ) can be expanded in powers of B. Apart from 
the different form of R^, this reproduces the Euler-McLaurin expansion. Expanding 
up to order B'^ we find, for instance, 

d ( W2\ .W2 127 /-p «N 



with Tq'^ = 2. This result is the counterpart of eq. ( p.lO|) which had been obtained 
with the exponential cutoff for which rg'^ = 1/6. In accordance with ( |2.15| ) we find 



an additional factor of 12 in the second term on the r.h.s. of ( |E.6| ). Hence also the 
value of the fixed point W2*{k) is 12 times larger than the result (|3.11|). In view of 



the discussion following ( |3.13D this means that for the mass-type cutoff Rk = Zkk^ 
this higher-order correction is much larger than for the exponentially decreasing R^ 
of ( |1.2|) . This is probably closely related to the ultraviolet problems and indicates 
that, though computationally more difficult to handle, the exponential cutoff (|1.2|) 
should be used for reliable estimates. 



Appendix F 

In this appendix we derive the flow equation in the truncation. We start from 
the evolution equation for ■w{i9) which follows from differentiating ( |3.5| ) with respect 
to 

d 

—w = (4 + r])id + Ai^w'^^^ 
at 

{oo 
m=l 

(2i?w2)m-i ^(^^g _ _ ^-^^^ ^ ^4 _ 2d)m -d+ y )?9u'^ + (4m - d)^ww''^^^ 
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+2(m - l){2^w^)"'~^w\w + 2^w){2mw + (4m - d)^w) 

4y;(3) + 279«;(4) 4:{w + dw^^^){3w + 2^w^^^) + 2^w{5w^^^ + 2^w^^^] 



(w + 2^wy 



{w + 2'dwY 
+ Mww^^^ 3 , 



§11] 



V,3 



w(w + 2'(9ti)) 2 it;2(w + 2'(9ti)) 



2rf 



+_rf(rf_l)(rf_2)(iV2-l)r2^ 



(3) d + 2 w 



w{w + 2'd'wy 

2' 



w 2 w"^ 

d,m/r,M\m—2 



+8{d - 2)vdg^ J2 - l)TmEmro'"'{2^) 



[F.l) 



771=1 



Evaluating this equation for 6 = yields eq. ( |3.18| ). Taking one further derivative 
at -i? = we obtain for d = 4 and k > knp the flow equation for ^3 (with C| — £'2 = 



29 ^4 7-1 



137 



10080 

d_ 
di 



and 



31 



30240 ^ ■ 

W3 = (8 + ri)w3 



+ 



9' 



167r2 



442 137 



315 210 



V ] T"3'^o'^ + ^(2 - r])T2r^"'w2 



i7 



4,2 



-(2 - r])r^{3wi - 5IW2W3 + 105w^) 
-3(2 - r/)(A^2 _ i)^4(^^^ _ Q^^^^ ^ ^2w^ 



(F.2) 



For the infrared cutoff (|1.2| ) one has = |,ro'^ 



— ^ and we observe that the 
perturbatively leading term ~ g"^ is negative. This implies a positive perturbative 
fixpoint value for the ratio w^/g'^. 

This discussion can easily be generalized for arbitrary Wn- In lowest order in g"^ 
the flow equation ( [bM] ) simplifies considerably {k > knp) 



d_ 
dt 



w 



Aw + Ai!}w^^^ 



m-2 



m=2 



This implies for n > 2 the flow equations 

^^wn = 4(n - l)wn -{d- 2)vdg'2^+'n\Tn4'' (c' - 
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(F.3) 



(F.4) 



and the infrared fixed point values 



in \ on— 1 

" ^ -nWy,'-{c'^-2E^) (F.5) 



id-2)vd9y^ 



n 



For d = 4 and with (|2.16|) this yields 



167r2 

(n^Ty(2n -1) V 2n 



dn = 7 TTTTT TTT^-^2n-2 B^n-l] (F.6) 



where for SU{2) 

do 127 rf, 221 



(F.7) 



T2 540' rg 37800 
The series of dn is alternating as long as the bracket is dominated by -1. 

The infrared fixed point in w^/ g"^ implies that the coefficient in ( |3.1|) diverges 

W2 ~ (F.8) 

and similar for W^. Because of the infrared divergence for /c — 0, a result of this type 
could never have been found in standard perturbation theory. It could, however, 
be derived using the "-^ ^"-approximation of the effective average action which 
we displayed in eq. ( |A.29D . For this purpose one can neglect the W^! terms on the 
r.h.s. of (|A.29|) and approximate Wl. = Wi. Then, with ( p.2|) inserted into ([A.29|) , 
one obtains for the /c-dependent terms in r^: 



Tk[A,A] = ^Tr,,iln 
-Tr,,. In 



WiVt + {wi - —)D ®D + RkiwiVr + {wi - —)D O D) 

O^k Oik 

-D^ + Rki-D^)] (F.9) 



If one extracts the F^^^-term from these traces, one finds a (renormalized) coefficient 
which equals exactly W2, as extracted from the fixed point ( p. 111 ). Clearly ( |3.8D goes 
beyond the "-^ ~ ^" approximation. The terms proportional to r] and to could 
not have been obtained in this approximation. 
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Finally, for k < knp the flow equation for follows from ( p^'.lD as 

d f 43/137 . 31 2\ 174 42 , 

di''' = + 16^ r^^°' ll05^^ + 315^0 + ^^^^ 

-6rt ( ^ - 17^ + 3541 

\Wi wf WiJ 
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